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Abstract 

We report measurements and searches for resonant — > K^h — ► K^jj de- 
cays where h is a 77, 77', r) c , ry c (25), XcOi Xc2, J/ip meson or the X(3872) particle. The 
results are based on a data sample containing 535 million BB pairs collected 
with the Belle detector at the KEKB e + e~ asymmetric-energy collider operat- 
ing at the T(4S) resonance. Signals are observed in the modes with 77 and 77', 
and we obtain evidence for a signal in the mode with 77 c . We measure B{B^ — ► 
K ± v _> k±7 7 ) = (0.87i^(stat)l°;$(syst)) x 1(T 6 , B(J3± - K+rj' - #±77) = 
(1.40+°;^(stat)+°;^(syst)) x 10" 6 and B(B ± -> i^ ± 7 ?c -> ^±77) = 
(0.22lg$(stat)±g;{£(syst)) x 10~ 6 . We set upper limits on the branching fractions 
of the other modes. 



Key words: B, X(3872), charmonia 
PACS: 13.25.Hw, 14.40.-n 
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1 Introduction 



We report searches for resonant B — > K ± h — > K ±r y r y decays, where h 
can be one of the following mesons: r], rj', r) c , r) c (2S), Xco, Xc2, J/ip or the 
X(3872) [1,2,3,4] particle. 

The nature and quantum numbers of the X(3872) particle are still being 
debated; based on analyses of the dipion mass spectrum [5,6] and angular 
distributions [5,7] for X(3872) -> n+n-J/ip, J PC = 1++ and 2" + are al- 
lowed. The 1 ++ assignment is also supported by signals observed in B — > 
(D°D 7i°)K [6] and in B -»• (D*°D°)K [8] under the assumption that they are 
indeed due to the X(3872) particle. The observation of X(3872) — > J/ipp [9] 
and X(3872) — > J/ip'j [10,11] indicates that C = +1. Evidence of a signal 
in the B^ — > K ± X(3872) — > if 77 channel would rule out J = 1 since the 
decay of a spin 1 particle (here the X(3872)) into two photons is forbidden by 
gauge invariance and Bose-Einstein statistics [12]. 

Many of the B^ — > K ± h and h — > 77 branching fractions involved in these 
decay chains have been already measured, as shown in Table 1. The B ± — > 
K r) and B^ K rj' modes are well established [14] and can be used as 
calibrations in the search for other B — > K ± h — > K ±, yy channels that have 
lower or unknown branching fractions. The B — > K ± J/ip channel can also 
serve as a control mode, since the J/ip is a spin 1 particle and cannot decay 
into two photons. 

The interference of B^ —>■ K ± rj c — > K ±r yy or B^ — > K Xco — * K ±/ ~fj with 
the radiative decay chain B — > K* ±r y — > K ±/ ~f / ~f can be used to measure the 
photon polarization in the b — > 57 quark transition [15]. Such measurement 
would provide a test of the Standard Model, which predicts the photon to 
be predominantly left-handed in b — > 57 decays and right-handed in b — > S7 
decays. The observation of the B^ —* K ± t] c — > K ±r ) r ) or B^ — > K Xco —* 
K ±, yy decay chain is the first step in this search for new physics, which could be 
achieved with about 10 ab _1 of data (thus requiring a Super B factory [16,17]). 
The non-resonant decay B^ — > K ±r y^ is very rare, with a branching fraction 
estimated to be of order 10 -9 [18] with a large background over the whole m 77 
phase-space from the resonant B^ —>■ K 7 — > K ±r yj channel [15]. 

In this study, we use a data sample of 492 fb -1 containing 535 x 10 6 BB pairs 
that were collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ (3.5 on 8 GeV) collider [19] operating at the T(4S") resonance. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists 
of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), a barrel- like arrange- 
ment of time-of-flight scintillation counters (TOF), and an electromagnetic 
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Table 1 

Current status of the measured branching fractions or 90% confidence level upper 
limits for B — > K^h and h — ► 77 (all values are taken from Ref. [14], unless 
otherwise indicated). The values in the last column are the expectations computed 
as the products B(B ± — > K ± h)xB(h — > 77). The decay chain B ± — > K ± h —* K ±r yj 
has only been observed for h = T). 



h 


B{B ± ^K ± h) 




B{B ± ^ I^h-tK* 


77) 


n 


(2.6 ±0.6) x 10~ 6 


(39.39 ± 0.24)% 


(1.02 ±0.24) x 


10 


-6 


ri 


(69.7 ± 2.8) x 10~ 6 


(2.12 ±0.14)% 


(1.48 ±0.11) x 


10 


-6 


Vc 


(9.1 ± 1.3) x 10~ 4 


(2.7 ± 0.9) x 10 -4 


(0.25 ± 0.09) x 


10 


-6 


Vc(2S) 


(3.4 ± 1.8) x 10" 4 


seen 








XcO 


(1.40i°;*j) x 10- 4 


(2.76 ±0.33) x 10~ 4 


(0.39 ± 0.08) x 


10 


-7 


Xc2 


< 2.9 x 10" 5 


(2.58 ±0.19) x 10~ 4 


< 7.5 x 10" 


•9 






(10.07 ±0.35) x 10- 4 


< 9.3 x 10~ 5 [13] 


< 9.4 x 10" 


•8 




X(3872) 


seen [1] 











calorimeter comprised of CsI(Tl) crystals (ECL) located inside a supercon- 
ducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return 
located outside the coil is instrumented to detect K\ mesons and to identify 
muons (KLM). The detector is described in detail elsewhere [20]. Two inner 
detector configurations were used. A 2.0 cm beampipe and a 3-layer silicon 
vertex detector was used for the first sample of 152 x 10 6 BB pairs (SVD1), 
while a 1.5 cm beampipe, a 4- layer silicon detector and a small-cell inner drift 
chamber were used to record the remaining 383 x 10 6 BB pairs (SVD2 [21]). 



2 Event selection and background rejection 



Kaon candidates are selected from charged tracks with the requirement C = 
£k/{£k + Ar) > 0.6, where Ck (At) is the likelihood for a track to be a 
kaon (pion) based on the response of the ACC and on measurements from the 
CDC and TOF. The kaon identification efficiency is between 84% and 90% 
depending on the ^77 signal mode with 7%— 11% of pions misidentified as 
kaons. Photon pairs are selected by requiring their energies in the laboratory 
frame to be greater than 100 MeV and their energy asymmetry A 77 = | ^ 7 |~^ | 

to be less than 0.9. We reject photons from 7r° decays by removing photon pairs 
with an invariant mass between 117.8 MeV/c 2 and 150.2 MeV/c 2 (2.5 standard 
deviations around the 7r° mass). We require a shower shape consistent with 
that of a photon: for each cluster, the ratio of the energy deposited in the 
array of the central 3x3 calorimeter cells to that of 5 x 5 cells is computed. 
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Table 2 



Nominal mass [GeV/c 2 ] 
mass windows [GeV/c 2 ] 


of the reconstructed particles and definition of invariant 
for photon pairs. 


Particle 


Mass 


Wide m 77 window 


Tight m 77 window 


V 


0.548 


0.4-0.7 


0.50-0.57 


>/ 


0.958 


0.8-1.1 


0.90-0.98 




2.980 


2.5-3.2 


2.82-3.05 


Vc(2S) 


3.637 


3.2-3.8 


3.44-3.70 


XcO 


3.415 


3.0-3.5 


3.25-3.50 


Xc2 


3.556 


3.0-3.8 


3.40-3.62 


J/ll> 


3.097 


2.5-3.2 


2.92-3.15 


X(3872) 


3.872 


3.0-4.1 


3.72-3.95 



The cluster associated with the most energetic photon of the candidate pair 
is required to have a ratio greater than 0.95 while the cluster from the other 
photon must have a ratio greater than 0.90 for the i? ± — > K ±j r) and JB ± — > 
K^rj' channels and 0.95 for the other channels. 

Pairs of photons are retained and associated to the corresponding meson when 
their invariant mass (?n 77 ) is inside one of the wide mass windows defined in 
Table 2. A mass-constrained fit of the photon momenta is performed to match 
the nominal [14] masses with the constraint that the photons originate from 
the interaction point. 

Charged B meson candidates are reconstructed starting from a kaon and a pair 
of photons, and they are selected by means of the beam-energy constrained 
mass, defined as M^ c = J -E^Lm — Pb an< ^ the energy difference AE = E* B — 
-^beam- I n these definitions, E^ cam is the beam energy and p* B and E* B are the 
momentum and the energy of the B meson, all variables being evaluated in 
the center-of-mass (CM) frame. We select 5-meson candidates with M\ )C > 
5.2 GeV/c 2 and -0.3 GeV < AE < 0.3 GeV. If more than one B candidate 
is reconstructed in an event, the best candidate is chosen by selecting the 
photon pair with the smallest x 2 of the mass fit, and if multiple kaons can be 
associated with this photon pair, the kaon with the highest C is chosen. 

The main background in all modes is due to continuum events, i.e. events com- 
ing from light-quark pair production (uu, dd, ss and cc). The rejection of the 
continuum is studied and optimized using a Monte Carlo (MC) sample having 
about 1.5 times the size of the data sample. Four variables are used to separate 
signal from continuum background: a Fisher discriminant based on modified 
Fox- Wolfram moments [22] , the B production angle with respect to the beam 
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Table 3 

Definition of the AE signal windows [GeV]. The M^ c signal windows are defined as 
Mbc > 5.27 GeV/c 2 for all modes. 



Particle 


AE window 


Particle 


AE window 


V 


-0.15 < AE < 0.10 


XcO 


-0.10 < AE < 0.10 


>f 


-0.15 < AE < 0.10 


Xc2 


-0.06 < AE < 0.06 


Vc 


-0.10 < AE < 0.10 




-0.09 < AE < 0.09 


Vc(2S) 


-0.08 < AE < 0.06 


X(3872) 


-0.09 < AE < 0.09 



in the CM frame, cos 6*, the flight length difference along the beam axis be- 
tween the two B mesons, and the flavor tagging information [23]. The Fisher 
discriminant, the B production angle and the flight length difference are com- 
bined into a likelihood ratio LR = L s j (L s + L u d sc ), where L s and L u d sc are the 
product of probability density functions (PDFs) of these variables for signal 
and continuum events. We use different LR cuts depending on the flavor tag- 
ging information. The continuum rejection is achieved by simultaneously opti- 
mizing the LR and m 77 cuts {tight m 77 window in Table 2) in order to maxi- 
mize the figure of merit in the signal windows (M^ c > 5.27 GeV/c 2 and AE as 
described in Table 3). The figure of merit is defined as S = N s /y/N s + N u d sc 
for the B ± — > K ± t] and B^ — > K^rf modes and e/y/ N u d sc for all the other 
modes, where N s and N UC [ SC are the expected number of signal and continuum 
events and e is the signal efficiency The expected numbers of events are com- 
puted for an integrated luminosity of 492 fb _1 and assuming the measured 
branching fractions [14]. 

Exclusive backgrounds from charmless B decays are studied using large MC 
samples having about 36 times the size of the data sample. In the B^ — > K ± r) 
channel, 56% of this type of background is from B — > K*rj with the rest 
being composed of several small contributions, the largest ones being due to 
B — > X s '~f and B^ — > ?77r ± . In the B — > K rj' channel, the dominant source 
(about 2/3) is B Xg'y, about half of which is from B — > K*{Ktt)^. For 
the other modes, about 95% of the charmless B decay contributions is due 
to B — > Xg'-y. The final state with K ± ir 0/ -f is a significant background for 
modes with charmonia and with the X(3872) resonance. It is suppressed by 
the requirement m Kl2 > 1.5 GeV/c 2 , where m Kl2 is the invariant mass of the 
system formed by the kaon and the lowest energy photon (in the laboratory 
frame) forming the K'j'-f candidate. For the B^ — > K rj c channel, the B — > 
-fT*(892)?7 c (77) background is the most relevant contribution. 

Another source of background is produced by the overlap of a hadronic event 
with a previous QED interaction (mainly Bhabha scattering) that has left 
energy deposits in the calorimeter. This off-time background is removed by 
using the timing information of the calorimeter clusters corresponding to each 
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Table 4 

Signal efficiencies for the two configurations of the detector. 



Particle 


e(SVDl) [%] 


e(SVD2) [%] 


Particle 


e(SVDl) [%] 


e(SVD2) [%] 


V 


15.8 ±0.1 


16.6 ±0.1 


XcO 


11.0 ±0.1 


11.6 ±0.1 


>/ 


14.6 ±0.1 


15.7 ±0.1 


Xc2 


10.4 ±0.1 


11.3 ±0.1 


Vc 


10.0 ±0.1 


10.9 ±0.1 




9.4 ±0.1 


9.7 ±0.1 


Vc(2S) 


10.9 ±0.1 


11.4 ±0.1 


X(3872) 


13.7 ±0.1 


15.0 ±0.1 



photon candidate. This timing information is only available for the most recent 
data, containing 258 x 10 6 BB pairs. For the rest of the data, we include 
the background in the fit described in the following section, by modeling it 
according to the off-time background events rejected from the most recent 
data. 



The tight m 77 windows overlap for some of the h decays, e.g. the mass window 
for rj c includes some candidates for J ftp and vice versa. Dedicated studies 
have shown that, for the dataset considered in this analysis, the only non- 
negligible cross-feed is due to i? ± — > K ± rj c events that are reconstructed in 
the B ± — ► K ± J/ip mode. This effect is included in the fit as described below. 



3 Fitting procedure and results 

We perform a two-dimensional unbinned extended maximum likelihood fit to 
Mb c and AE. The signals are described using PDFs modeled with the product 
of a Crystal Ball function [24] for M^ c and three Gaussian functions for AE, 
while the continuum background is modeled with an ARGUS function [25] for 
Mbc and a first order polynomial function for AE. The effect of neglecting the 
correlation between Mb c and AE has been studied using MC signal samples 
embedded in toy continuum samples; the number of signal events returned by 
the fit is found to be 1-3% smaller than the true number, depending on the 
h mode. We take this bias into account by correcting the signal efficiencies 
and adding a systematic uncertainty. Table 4 lists the corrected efficiencies 
obtained for each mode in the two sub-samples with different inner detec- 
tor configurations. The signal PDF parameters are determined on MC signal 
events. The M^ c resolution and the AE resolution and mean are then cor- 
rected using a control sample of B ± — > K ± n° events. The BB and off-time 
backgrounds are modeled with two-dimensional KEYS [26] PDFs extracted 
from MC events and from the off-time data sample, respectively. The nor- 
malizations of the BB and off-time backgrounds are fixed in the fit. For the 
B ± -> K^J/ij) mode, the K±i) c cross-feed is included with normalization fixed 
to the value obtained in the corresponding signal fit. 



S 



Table 5 

Systematic uncertainties on the signal reconstruction efficiency. 



Source 

t. J ' > 111 V . V 


Uncertaintv ynl 


Pnntnn vprrin t;t vi 1 cf ion PTTir'ipnr'ir 


2 x 2.2 


Tracking efficiency 


lxl 


Kaon identification efficiency 


2 


m 77 cut efficiency 


3.6 


LR cut efficiency 


6.9 


MC statistics 


1.0 


Fit bias 


0.5 


Total 


9.3 



The fit is performed for M^ c greater than 5.2 GeV/c 2 and for AE between 
—0.3 GeV and 0.3 GeV. The likelihood is defined as: 

£ = x n(E^(Kc AS*)) (1) 

where i runs over all events, j runs over the possible event categories (signal, 
continuum background and other backgrounds), Nj is the number of events 
in each category and Pj is the corresponding PDF. 

The data are divided into sub-samples based on the SVD configuration and 
the availability of the timing information needed for the rejection of off-time 
background. 

The fit variables are the branching fraction (£>) and the continuum background 
normalization and PDF parameters, except the ARGUS endpoint which is 
fixed to -Ebcam = 5-289 GeV. The number of signal events is then defined as 
S k = B x e k x N^y- where Nt-s is the number of BB events and e k is the signal 

tit) JDJD 

efficiency, both evaluated for sub-sample k. 

The branching fraction obtained from the fit depends on the following param- 
eters that can give rise to systematic uncertainties: 

(1) parameters related to particle reconstruction and identification and to 
signal selection, which affect the signal in a very similar way for all h 
modes, as summarized in Table 5, 

(2) signal PDF parameters (0-5% uncertainty), 

(3) normalization of the charmless B and off-time backgrounds and of the 
K ± r\ c cross-feed for the B ± —> K ± J/ijj mode (1-10%), 

(4) number of BB events (1.3%). 
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Systematic uncertainties related to the m 77 and LR requirements are evalu- 
ated by comparing efficiencies in data and MC using a B — > K 7T° control 
sample. Systematic uncertainties are included in the likelihood function by 
integration. The statistical likelihood is convolved with the probability distri- 
bution of the systematics parameters listed above, computed as the product of 
Gaussian terms, one for each parameter. A MC integration is performed over 
the phase space of the systematics parameters, yielding a new likelihood func- 
tion, £ syst , that includes all systematic uncertainties. The fit results quoted 
below are all extracted from £ syst . The central value B is the B at which £ syst 
has its maximum and the errors 8^ t are defined by : 

Jr„_i_A- ^syst «£> 
Bo+S tot =06g (2) 



Jo ^syst dB 



where the integration interval is chosen such that all points outside the interval 
have a lower likelihood than those inside. The positive (negative) systematic 

error is computed as iy^t — ^stat where 5jJ at is the positive (negative) 
statistical error. The significance of the measurement of the branching fraction 



is defined as J 2 (In C syst (B = Bo) — ln£ syst (£> = 0)). For modes in which no 
significant signal is found, the 90% credible upper limit, Sumit, is computed 
using a Bayesian approach with a flat prior, according to: 

j ° 1 ^ syst m = 0.9 (3) 

Jo ^syst dB 



The fit results for all modes are summarized in Table 6. We observe signals in 
the B ± — > K ± r] and B^ 1 — > K ±q r( modes and obtain evidence for a signal in the 
B^ — > K ± rj c channel, while we see no signal in the other modes. We report the 
first measurements of B^ — >• K ± rj' and B^ —* K ± rj c channels in the K ±r yy final 
state. We measure B(B ± -> -> ^77) = (0.87lg;}f(stat)lg;S?(syst)) x 

10 -6 in agreement with Belle's measurement of this mode with the same 
dataset [27], B(B ± -> K^rj' -> ^77) = (1.40t[J;^(stat)±^(syst)) x 10" 6 
and B(B ± -> -> ^±77) = (0.22lg$(stat)l^(syst)) x 10" 6 . All mea- 

sured branching fractions agree with the values shown in the third column of 
Table 1. Fit projections are shown in Figures 1 and 2; in each plot the variable 
that is not shown is restricted to be in the signal window. 

For the modes where no significant signal is observed, we extract the following 
90% probability upper limits: B(B ± -> K^J&S) -> ^77) < 0.18 x 10~ 6 , 
B(B ± -> K±Xda ^ ± 77) < 0.11 xlO" 6 , B(B ± -> K ± Xc 2 -> ^77) < 0.09 x 
10~ 6 , B(B ± -> i^J/V -> ^77) < 0.16 x 10" 6 and B(B ± -> AT ± X(3872) -> 
^77) < 0.24 x 10~ 6 . Whenever the branching fraction of B^ — > A' =l= /i has 
been measured elsewhere, we also perform the fit by constraining B(B ± — > 
K ± h) to the measured value [14], thus extracting an upper limit on B(h — > 77). 
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Table 6 

Signal yields, branching fractions and significances (5) results for — ► K^h — ► 
i^ ± 77. The first uncertainty is statistical, the second one is systematic. Limits are 
calculated at 90% probability. 



Resonance Yield B(B ± -> K±h -» #±77) (10~ 6 ) 5 £(7t -> 77) 



1/ 

%(25) 


7fi +14 

114 ±13 

13.3+tf 
4 0+ 39 


n o 7 +0.16+0.10 
u -°' -0.15-0.07 

1 An +0.16+0.15 
i -^ u -0.15-0.12 

n 99+0.09+0.04 
u - zz -0.07-0.02 

< 0.18 


7.3 
13.8 
4.1 


(0 nl +0.23+0.23\p7 
V z - ul -0.22-0.19J /l 

(2.4^|1°;I) x 10- 
< 8.2 x 10~ 4 


XcO 




< 0.11 




< 9.5 x 10~ 4 


Xc2 


-0 3+ 2 6 


< 0.09 






J/lj> 


4 +2.8 


< 0.16 




< 1.6 x 10" 4 


X(3872) 


-0.9^-2 


< 0.24 




< 1.1% 



The uncertainty on B(B ± — > K h) is included as a source of systematic 
uncertainty. We obtain B(x c o — > 77) < 9.5 x 10~ 4 , B(r] c (2S) — > 77) < 8.2 x 
10~ 4 and B(J/ip — > 77) < 1.6 x 10~ 4 at 90% probability. Similarly, for the 
B ± —> -R' ± ?7 C mode, we determine B{r\ c — > 77) = (2.4±o|(stat)io I(syst)) x 
10" 4 . 



The absolute branching fraction B(B ± — > K ± X(3872)) has not yet been mea- 
sured. However, there are measurements of the product of this quantity and the 
branching fractions of different decays of the X(3872). Assuming that X(3872) 
decays to J/?/>7r + 7r~, J / ilm + ix~ ix° and J/ipj saturate all possible decays of the 
X(3872) and taking the values of the corresponding products from [14,10,11], 
we derive a conservative upper limit £>(X(3872) — > 77) < 1.1% at 90% prob- 
ability. 



4 Conclusions 

A search for resonant B ± — > K ± h — > K ±r yy decays, where the resonance 
h can be 77, 77', r/ c , r) c (2S), Xco, Xc2, J/i> or X(3872), has been performed in a 
sample containing 535 million BB pairs. We have observed B^ — > K ± rj and 
B^ — * K ± r]' with significances of 7.3 and 13.8, respectively, and we have 
obtained evidence for 5 ± — > K ± rj c with a significance of 4.1. No evidence of a 
signal is observed in any of the other modes and 90% probability upper limits 
are set on the corresponding branching fractions. The measured branching 
fraction for B ± — > K ± rj — > K ± ^^ is in agreement with Belle's measurement 
of this mode with the same dataset [27]. We report the first observation of 
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— > K^rj and the first evidence of B ± — > K rj c in the K ±r ) r ) final state. 
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Fig. 1. Mb c and AE projections together with fit results. The first row presents 
the i? ± — ► K^f] mode, the second one Z? ± — ► K^rj ', the third one — ► K^r\ c 
and the last one — ► K^r] c {2S). The points with error bars represent data, the 
thick solid curves are the fit functions, the thin solid curve is the signal function, 
the dashed curves show the continuum contribution and the dotted curves show 
the sum of all background contributions. The hatched area present in the whole 
AE region is the contribution from the charmless B decays. The hatched area 
around AE = —0.2 GeV in B^ — > K^rj (-B ± — ► K^ri c ) shows the contribution 
from B —> K*7] decays (B —* K*rj c ). The filled area around AE = 0.05 GeV 
in the B^ — > K^rj plot is the contribution from B^ — ► ir^r/. The filled area in 
B^ — ► K^rj c {2S) is the contribution from the off-time background. 
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Fig. 2. Mb c and AE projections together with fit results. The first row presents the 
i? ± — ► K^XcO mode, the second one i? ± — ► K ± \c2, the third one — > K^J/i/i and 
the last one B ± -> K ± X(3872). The points with error bars represent data, the thick 
solid curves are the fit functions, the thin solid curve is the signal function, the dotted 
curves show the sum of all background contributions, the dashed curves show the 
continuum contribution, the hatched areas are the contribution from the charmless 
B decays and the filled areas the contribution from the off-time background. In the 
B ± -> K+J/ip plots, the B ± -> K^c cross-feed is visible in the thin solid curves 
as a small peaking background in M]~> c that is concentrated around 120 MeV in AE. 
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